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>. ■ Abstract 

Using a recently proposed measure for divisibility of a dynamical map, we study the non- 
Markovian character of a quantum evolution of a driven spin-S 1 system weakly coupled to a bosonic 

od 

bath. The complete tomographic knowledge about the dynamics of the open system is obtained by 
the time-convolutionless master equation in the secular approximation. The derived equation can 
be applied to a wide range of spin-boson models with the Hermitian or non-Hermitian coupling 

X . 

operator in the system-environment interaction Hamiltonian. Besides the influence of the environ- 
ed \ 

mental spectral densities, the tunneling energy of the system Hamiltonian can affect the measure 
of quantum non-Markovianity. It is found that the non-Markovian feature of a dynamical map of 
a high-dimension spin system is noticeable in contrast to that of a low-dimension spin system. 
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I. INTRODUCTION 



The realistic dynamics of an open quantum system is expected to deviate partly from the 
Markovian evolution which is idealized under the influence of memoryless environments IrtZ]. 
The coupling of a system to an environment inevitably leads to the loss of information and 
dissipation of energy js-lO]. In the Born-Markov approximation, the decoherence process 
shows a monotonic decrease with time. However, many kinds of condensed matter systems 
will undergo non-Markovian evolutions because of complicate interactions with their sur- 



roundings llNl3|. The non-Marokovian dynamics of an open quantum system can play an 
important role in a variety of modern applications such as quantum optics , solid 



state controllable technology 



15], quantum chemistry 



16| . and quantum information jl7 |. 



Under the influence of environments with memory effects, the non-Markovian dynamics may 
preserve the existence of quantum entanglement for longer time 18|, ll9[. Therefore, it is 
necessary to efficiently quantify the measure of non-Markovianity. 



Recently, several approaches have been put forward 2014231] . The measure based on the 
distinction of quantum states is used to describe a temporary flow of information from the 
environment back to the system. The non-Markovian behavior will be remarkable if the trace 



distances of evolving states increase 



20|. This measure emphasizes the physical features of 



the system-environment model. According to the composition law, another measure was 
proposed to quantify the deviation from the divisible map. An indivisible map shows the 
emergence of the non-Markovian dynamics. The mathematical description of this measure 
is simple due to no optimization procedure 2jJ, |22|. Meanwhile, some operational measures 
were offered, including the measure by assessing the increase of entanglement between an 
open system and an isolated ancilla j2l| and the measure by the quantum Fisher information 



flow 



23] 



Many theoretical models of efficient two-level systems coupled to their environments 



were investigated to quantify the measure of non-Markovianity of the evolutions 



24-3(J. 



As is known, hig 



in the nature 



i-dimension spin systems or effective <i-level systems extensively consist 



31 



33| . It is interesting to understand the complicate evolutions of these 



open systems. On one hand, we perform a reasonable unitary transformation to establish 
the second-order time-convolutionless master equation for a driven spin-S" system. With 
the consideration of the interaction Hamiltonian between the system and environment, the 



master equation is valid for both Hermitian and non-Hermitian system coupling operator 
cases. This is one point of the motivations of our present work. On the other hand, we utilize 
the measure for the divisibility to quantify the the non-Markovianity of the dynamical map. 
Through the strict theoretical analysis, we study in detail how the dimension of open system 
affects the degree of quantum non-Markovianity. This is the other point of the motivations 
of our investigations. 

The article is organized as follows: In Sec. |TT1 we briefly describe a model of a driven spin- 
s' system coupled to a bosonic bath. In the eigenrepresentation of the spin-S Hamiltonian, 
the evolution of the open system is obtained perturbatively up to the second order in the 
weak interaction between the system and bath. The divisibility of the dynamical map is 
introduced and calculated to quantify the degree of the deviation from a Markovian process 
in Sec. IIHI Finally, in Sec. IIV} we summarize our results and draw some conclusions. 

II. EVOLUTION EQUATIONS OF A DRIVEN SPIN-S SYSTEM COUPLED TO 
THERMAL BATH 

A. Hamiltonian of system-environment model 

Our system-environment model consists of an arbitrary driven spin-S system which is 
weakly interacting with a thermal bosonic environment. The general total Hamiltonian is 
composed of three parts, 

H = H S + H E + Hj , (1) 

where Hs = cu s n • S describes the Hamiltonian of a driven spin-S* system. The parameter 
u s denotes the transition frequency between any two neighboring energy states. The unit 
vector n = (sin a cos sin a sin (p, cos a) is used to characterize the tunneling element of 
the Hamiltonian. These angles satisfy that a G [0, tt) and (p G [0, 2tt). The operator of 
spin-S is expressed by three components of S = (S x , S y , S z ). Using the raising and lowering 
operators S± = S x ± iS y , the Hamiltonian of a driven system is also written as, 

H s = -u) s [mna(§ + e-* fi + S_e iLp ) + 2cosaS,] . (2) 

S±\m) = a/(S ± m + 1)(S =F m)|m±l) where \m) is the eigenvector of S z with the eigenvalue 
m. Here, the tunneling energy is mainly determined by the angle parameter a. 
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To conveniently investigate the dynamics of a driven spin system, we turn H$ into the 
form in the eigenrepresentation spanned by the basis {\m),m = S, . . . , — S}, 



H s = U ] HsU 



(3) 



where the unitary transformation is expressed as, 

s 



U 



m=—S 



( „(S) .(s-i) 



a 



s 



s 



(S) (S-1) 
a S-l a S-l 



a {s) a (s " 1} 

a -S a -S 



7 (- s ) \ 

l S-l 



(-s) 



(4) 



Here the normalization state \ip m ) = Ylk a ^\^) ^ s ^ ne eigenvector of the Hamiltonian H$- 

In regard to weak couplings of an open system to an environment, we may reasonably 
consider the environment as a collection of harmonic oscillators. The Hamiltonian of the 
heat bath is written as, 

//, ( 5 ) 

j 

where fot and bj are, respectively, the creation and annihilation operator of the jth harmonic 
oscillator with the frequency Uj. 

The interaction Hamiltonian between the system and environment is expressed as, 



(6) 



The two cases of the system coupling operator ft are considered here. In one case, the 
operator is Hermitian, i.e., IF = EL The other case of a non-Hermitian operator El^ ^ II 
is also involved. The weak coupling strength gj can be characterized by the environmental 
spectral density J(u>) = J2j Ifl'il 2 ^^ — u j) M- 



B. Quantum evolution equation 



On the assumption of the weak system-environment couplings, the time-convolutionless 
master equation for the second-order perturbation after the unitary transformation U is 
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written in the interaction picture to describe the dynamics of the open system, 

Jt m = ~^J dt ^E[Hi(t)AHi(ti),p(t)®PE)] , (7) 



where the density matrix of the open system p(t) = Tr^[p tot (t)] is obtained by partially 
tracing over the freedom degrees of the environment. The initial total state is assumed 
to be factorized as ptot(O) = p(0) ® Pe where pe = exp(— H E /ksT)/ Tr[exp(— H E /k B T)} 
denotes the thermal equilibrium state of the environment. Here the Planck constant h and 
Boltzmann constant k B are taken to be one. 

The general interaction Hamiltonian in the interaction picture is expressed as, 

Hj(t) = U^e^H^-^U 

= E E ^d-e^-^^rnb] + h.c. , (8) 

j l,m 

where the Hamiltonian H represents H = H$ + H B - The transition operator is defined as 
a l,m — 10 (fn\ and the notation h.c. denotes the Hermitian conjugate of the former term in 
Eq. ([8]). The coefficients c/ >m are calculated by U^IUJ = J2i m c i,m,o'i,m- 

In the field of condensed matter physics, many kinds of effective d-level quantum sys- 
tems or some high- dimension spin systems have their coherence time tq which may be much 
smaller than the correlation time of the thermal environment te [31| . In the condition 
of Ts <C te, we may reasonably adopt the secular approximation by neglecting the high- 
frequency oscillating terms. This secular approximation is almost equivalent to the rotating 
wave approximation in quantum optics 8j. The necessary derivation of the quantum evo- 
lution equation is shown in the appendix. The general expression of the quantum master 
equation for the dynamics of the driven spin-S 1 system is presented as, 



-~ P {t)= 53(r w + f, 
i 

+ E( A ^ + A^)c fc>fc c^ ■ a-i t ipa k ,k + h.c. 



( ]fi <-' / - '■' ~i,i)\ci,m\ • \?i,ipvi,i — 2 a i,iP ~ 2 f}CTl '^ 
I 



l^k 

+ ^2(Xl,m\Cl,m\ 2 + Tm,l\ C m,l\ 2 ) ' (p m l,P°l,m ~ \°l,lP ~ \p°l,l) ■ (9) 
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The decay rates in Eq. (Q are obtained as, 



A/,m(*) = J dh J dco J(w)re <[w - (, - m)w ' 1( *-* l) 
Ai, m (t) = f Q dt ±f ^^(w)(^+l)e ![w " (, " mK]Nl) , 



where the other decay rates T/ jm = 2Re(A/ >m ) and Y^ m = 2Re(A^ m 
Re(-) represents the real part of a complex number, r = l/(e W//T 
number of the thermal environment at a finite temperature. 

For the simplest example of S = |, the coefficients for the Hermitian coupling operator 
IT = S. are calculated by ci 1 = — c=i =1 = 1 cos a and Ci =i = c*_ x x = |e iv * sin a. Then, 

2 ' 2 2 ' 2 Z 2'2 ~2~'2 

the evolution equation just describes the dynamics of famous spin-boson models. In the 
other case of non-Hermitian operator TL = S_, the coefficients are also obtained as ci i = 



The symbol function 
1) signifies the mean 



— c=i =i 

2 ' 2 



^e t!p sin a, ci -i 

2 2 ' 2 



— e 2tip cos 2 % , and i 

2 2 ' 2 



— sin f. This quantum master 



equation can characterize the dynamics of the spin system coupled to the environment in 
the form of the dipole-dipole interaction. This result coincides with the quantum master 
equation in Ref. [28]. 

It is also interesting to present the quantum master equation for the dynamics of a driven 
spin-1 system. The unitary transformation for the case of S = 1 is given as, 



U 



( 



\ 



cos 2 § e-*? 



sin ae 



•tip 



sin 2 f e"^ 



sin a 



cos a 



sin a 



sin 2 f 



v/2 



sin ae 



up 



- cos 2 f 



/ 



(10) 



For the example of II = S z , the coefficients in Eq. (jUJ) are written as cy = I ■ cos a, q j0 = 
= — £ • sin a, (/ = 1, —1), and c\ _i = c^i = c 0; o = 0. In this case, the expression of 
the quantum master equation is simplified as, 



d „, . 



7o ■ - a-\,-i)p{(Tx,i - o"-i,-i) - - - - °"-i,-i)] 



1 



1 



+ 7+ ■ [cTi i0 pcro,i + o-Q,-ipo--i,o - 2°"o,oP - 7^00,0 



1 



1 



-cr-i-ip- -p^-i-ij 



1 1 1 1 

The decay rates in Eq. (iTTj) are written as 70 = [ci 3 i | 2 (ri ; i + T^i), 7+ = |co,i | 2 (r ,i + ri,o), 
and 7_ = |ci i0 | 2 (r li0 + f 0)1 ). 
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III. NON-MARKOVIANITY MEASURED BY DIVISIBILITY 



We will quantify the quantum non-Markovianity based on the divisibility of a dynamical 



map introduced by Rivas, Huelga, and Plenio (RHP) 21] . According to the composition 



law, a trace-preserving completely positive map Q(t, 0) is divisible if it can be expressed as, 

e(t + h, o) = e(t + h, t)e(t, o) , (12) 

and Q(t + t%,t) is completely positive for any t , t\ > 0. The divisible map Q(t + £i,0) is 
Markovian exactly. It is shown that the quantity satisfies, 

9lt) = Inn II^»8(' + '..«WH.-1 > , (13 ) 

ti-»0+ ti 



where 1^ is the identity matrix in the ci-dimension Hilbert space T-i d and the density matrix 

D + 
d 



Pf signifies the maximally entangled state in the product Hilbert space 1-L d ® d . The function 



li represents the norm of the matrix A. The value of g(t) is strictly positive if and only 
if the dynamical map is indivisible. Therefore, the condition of g(t) > determines the 
emergence of the non-Markovian dynamics. The quantity g(t) gives rise to a measure for 
quantum non-Markovianity defined as, 

N RHP (&) = -—, 1= dtg(t). (14) 

-L + 1 Jo 

For the dynamical map of an open spin-S* system, the maximally entangled state is written 
as P^ = 2 ^2j n =-S ® "where d = 2S + 1. When the dynamics of the open 

system is characterized by the quantum master equation d t p = C t p, the key quantity g(t) is 
expressed as, 

u\ v IKMi)||i ~ 1 
g{t) = km , (15) 

tl-5>0+ t\ 



where the matrix e(t, tt) = [(I d + t x C t ) <g> IJP/ [21]. 

For the case of S — 1, applying the dynamical map like Eq. ( ITT]) , we can obtain the 
matrix e(t, t\) as, 
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1 

£= 3 



1 - *!7_ 
ti7+ 

£i7_ 

1 — tix/ 






(16) 



1-hy 







l-t 1 (7+ + 7_) 

*!7 + 



ti7_ 

1 — 

where the elements are defined as y = y + y +7_, t> = y + -y + 270, and w = y + y- + 7+- 
The key quantity is calculated as, 

9® = A E $ ^') > 
d i=o,± 

where the value of the function is given by = x, (x < 0) and = 0, (x > 0). 

To analytically study the effects of the dimension of the open system on the non- 
Markovianity, we can offer the general expression of the matrix e(t,ti) as, 



l-t x v \ 




1 — t\W 




1 - *i7+ / 



e(t,ti) = + %{^2(ri,i + ^i,i)\ci,i\ 2 (Sij5 n iai,i - -5ij(Ti n - «<Wj,«) 



3,n,l 



+ ^ ( A M + ^-l,l) c il c k,kSlj6 nk (Tl,k + h.c. 
j,n,l^k 



+ \ci,m\ 2 Ti tm (5ij5 n i(J mtm — ~5lj(Tl n — ~Snlffj,l) 



(18) 

where 5y = (5^ is a delta function. For the case of Hermitian coupling operator II — S z , the 
coefficients in Eq. ( fl"8l) are calculated as = Z ■ cos a, / = —5, ■ ■ • , S. With respect to the 
special case of a ~ 0, the decay rates T/ / = r and = T dominate in the dynamics of 
the open system because the off-diagonal elements \c^ m \ 2 ~ 0, (I ^ m) . In this condition, 
we can analytically obtain the quantity g(t) as, 

/(d) 



d 



cos 2 a[$(r ) + $(r )] 



(19) 



s 



where the quantity g(t) is closely related to the dimension of the system and the function 
f(d) satisfies f(d+3) = 3[f(d+2)-f(d+l)] + f(d) + l where /(2) = 1, /(3) = 4, /(4) = 10. 
It is straightly seen that the degree of non-Markovianity can be increased with the dimension 
of spin system. 

IV. CONCLUSIONS 

We derive in detail the second-order time-convolutionless master equation for the dy- 
namics of a driven spin-5 1 system weakly coupled to the thermal bosonic environment. The 
quantum evolution equation expressed in the eigenrepresentation of the Hamiltonian H$ is 
valid for both Hermitian and non-Hermitian coupling operator in the interaction Hamilto- 
nian Hj. We also present the quantum master equation for the case of S = 1 in the secular 
approximation. The quantum evolution equation can be generally applied to the dynamics 
of any effective (i-level quantum system. To conveniently characterize the degree of the de- 
viation from Markovian dynamics, we analytically study the indivisibility of the dynamical 
map of a driven spin-S* system. The case of S = 1 is also considered. The indivisibility of 
the map is mainly determined by the negative values of the decay rates and also related to 
the tunneling energy. Trough the theoretic analysis, we obtain the interesting result that the 
non-Markovian behavior of high-dimension systems may be remarkable compared to that of 
low-dimension systems. 
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Appendix A: Derivation of quantum master equation ([9]) 

In this appendix, we derive in detail the second-order two-convolutionless master equa- 
tions (J7J) and ([9]). Firstly, the interaction Hamiltonian Hi(t) can be rewritten as, 



Hj(t) = B\t) <g> C(t) + B(t) g> C\t) , 
9 



(Al) 



where the operators are defined as B = J2j9j e W]t ^j and C = Ylim c hm e ^ m ' Uat oi,m- The 
dynamics of the total state of the system-environment model is calculated as, 

7 • -| ft 

-p tot = -^[Hj(t),p tot } --J dhlHjit), [^j(ti), Aot^i)]] . (A2) 

On the assumption of the factorized initial total state p(0) ® Pe, we may safely make the 
approximation of Tr E[Hi(t)pE[ = in order to eliminate the first term in Eq. (1A2j) [3]. 
When the Born approximation is adopted, we can replace p to t(£i) with p(t) <8> pe to obtain 
Eq. ([7j) after performing the partial trace over the freedom degrees of the environment. 
Then, the second-order time-convolutionless master equation can be expanded as, 

d - p{t) = f dt x [ duJ{u)r ■ e^-^lC^hYpCit) - C(t)C\h)p] 



dt 

+ J dhj duJ{u){r + l)-e iu{t ~ tl \C{t)pC\t 1 )-pC\t 1 )C{t)] + \i.c. (A3) 
The term of C^(ti)pC(t) in Eq. ( )A3j) is obtained as, 



ct(fi)pc-(f) = EE^/ i|(l " m)tl+H)tK ^-i • 

i,m k,j 

When the secular approximation oik = l^j = m and k = j ^ I = m is made, we can 
simplify the above equation as, 

C^^pCit) = cliCk,m,ip(Th,k + hrnl^-^^^pa^ . (A4) 

l,k Ijtm 

Similarly, other terms in Eq. (|A3[) are written as, 



l,k l^m 
l,k l^m 

p&WC® = tikcliWVij + hml'e-^^-^pa^m , (A5) 

l,k l^m 

where 5u = dik is a delta function. By substituting Eqs. flA4j) and (IA5j) into Eq. ( 1A3j) . we 
can obtain the quantum master equation of an open spin-S 1 system like Eq.(|9]). 
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